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Abstract 

We consider a model of static cosmic string loops in type IIB string theory, where the 
strings wrap cycles within the internal space. The strings are not topologically stabilised, 
however the presence of a lifting potential traps the windings giving rise to kinky cycloops. 
We find that PBH formation occurs at early times in a small window, whilst at late times we 
observe the formation of dark matter relics in the scaling regime. This is in stark contrast to 
previous predictions based on field theoretic models. We also consider the PBH contribution 
to the mass density of the universe, and use the experimental data to impose bounds on the 
string theory parameters. 



1 Introduction 

In recent years there has been a renewed effort to test string theory in a cosmological context. 
This is due in part to the availability of increasingly precise data from experiments such as 
WMAP [1] and SDSS [2]. However, it is also due to theoretical advances which have resulted in 
a better understanding of the compactifications of the theory down to (3 + l)-dimensions. Since 
such compactifications are typically warped this means that mass scales in the effective theory 
can be significantly reduced. One important consequence of this is that superstrings may have 
a much smaller tension than first realised [31 [3] • 

Originally Witten [5] ruled out the notion that i^-strings could play the role of cosmic- 
strings because of their extremely high mass density. The un-warped string tension is so large 
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that the presence of cosmic F-strings (formed after inflation) would be immediately evident in 
the CMB. This effectively killed the field until GKP (Giddings-Kachru-Polchinski) [B] showed, 
in the context of type IIB strings, that by turning on non-trivial fluxes threading cycles in a 
class of compact manifolds one could obtain highly warped four-dimensional backgrounds. The 
warping acts in such a way as to reduce the overall tension of any object, and therefore it was 
possible to evade the observational bounds. 

Simultaneously there has also been renewed interest in models of open string inflation. In 
such models the energy density of the inflaton field (s) is provided by the geometric distance 
between the D3-brane of our universe and parallel branes or anti-branes. Since such branes/anti- 
branes are charged under massless RR-lorm fields it is expected that large numbers of F and 
D-strings, as well as bound states of multiple strings, will be produced at the end of inflation 
[H [71 [8]. Generically these objects are called (p, g')-strings, since they are typically a bound 
state consisting of p F-strings and q D-strings [9l \T0\ [TT] . Observationally however this raises 
a potential problem as these strings are not red-shifted away and therefore (presumably) fine 
tuning is required to restrict their number to 0(1) per Hubble radius at the present epocHj. In 
this paper we will simply assume that some phase of open string inflation has occurred, though 
it remains an open problem to generate such a configuration in an explicit model of string theory 
inflation. 

Given that these strings exist in a higher dimensional theory, one could also imagine that they 
wrap cycles within the internal space [12] . Such cycles could be either " smooth" or " lumpy" from 
a four dimensional perspective. A string that wraps a series of internal cycles at separate points 
in four-dimensional space would appear as a necklace. That is, as a system of monopoles or 
"beads" connected by string segments. Alternatively, strings which wrap the compact directions 
smoothly along their four dimensional length would appear to have a continuously varying mass 
density. The exact nature of the variation depends on the geometry of the compact space, 
though in general we would expect the effective tension to be periodic. 

Matsuda has proposed that a necklace structure may form from a smoothly varying config- 
uration as the string relaxes to a quasi-stable state [13]. In this case the periodic variation of 
effective string tension is viewed as the sum of the standard string tension plus a lifting potential 
in the (angular) compact directions. If the angular directions are flat the energy of the string is 
minimised for a smoothly varying configuration, resulting in a constant effective string tension, 
but in the presence of a potential a necklace structure is energetically favoured. Furthermore the 
existence, or absence, of a potential affects the stability of the extra-dimensional windings when 
strings chop off from the network to form loops. In the absence of a lifting potential windings 
must be stabilised topologically giving rise to objects called "cycloops". This is not true for 
necklace solutions which may exist, at least over cosmologically relevant time scales, even if the 
compact space is simply connected. The quasi-stability of the necklace solution is discussed in 

''This is necessary for the string embedding considered here. We argue for a UV completion of cosmic string 
scenarios so we must consider the origin of the strings/necklaces in the parent theory. Since the most well 
modelled inflationary scenarios in type IIB string theory are brane inflation scenarios, we will take this as the 
inflationary mechanism. Ifowever brane inflation typically results in the formation of large numbers of defects 
after the inflationary epoch. Since we do not observe such objects in the visible horizon we must allow for some 
kind of tuning of the theory prior to inflation which ensures that (in the post inflationary epoch) the relative 
numbers of branes, anti-branes and fluxes cancel leaving a small number density of residual defects. 
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greater detail in section 5.5. 

The cosmological consequences of cycloops were first investigated by Avgoustidis and Shel- 
lard [11]. They showed that, unhke an ordinary string loop, a decaying cycloop leaves a topo- 
logically trapped remnant when it reaches zero radius. This remnant appears as a monopole to 
a four dimensional observer and may be interpreted as a dark matter particle. Like ordinary 
string loops, cycloops also have a small probability (/ ~ 10"^*^) of collapsing to form PBH's in 
the course of their first oscillation. These PBH's may then decay to leave topologically stable 
Planck-mass relics, though this possibility has not been thoroughly investigated [T^ [T5]. 

By contrast Matsuda has claimed that a large fraction / ~ 1 of all necklace loops eventually 
collapse to form PBH's via a separate necklace-specific process. This claim is based on the 
assumption that mass-energy may only be lost, via gravitational wave emission, from the four- 
dimensional string segments, leaving the bead mass unchanged. Thus when the necklace loop 
reaches its minimum radius it will undergo collapse if the mass of the beads is large enough 
to produce a Schwarzschild radius greater than the string widthH He therefore proposed that 
small loops produced at very early times may form stable relics which also act as dark matter 
candidates. Conversely he argued that loops created at later times would be larger and hence 
likely to contain enough mass in their beads to cause them to collapse into PBH's |15j . 

However this original schematic analysis involved a number of simplifying assumptions, such 
as the existence of a time-independent lifting potential and hence a constant bead mass for 
necklace loops formed at different epochs. The initial inter-bead spacing was also assumed to 
scale like the entropy distance at the time of network formatiorl^l and the dynamical evolution 
of the inter-bead distance was modeled by the standard string-monopole network evolution 
equations, originally proposed by p!7] . 

In the following analysis we attempt to construct a more concrete model of necklace forma- 
tion based on ideas from type HB string theory. We calculate the explicit form of the lifting 
potential for a loop of string with extra-dimensional windings in the Klebanov-Strassler geom- 
etry \18\ I19j. Using realistic models of winding and loop formation we see that potential itself 
evolves dynamically resulting in a time dependent bead mass. This shows that the first of these 
assumptions must be modified, at least for certain backgrounds. 

Additionally the decay signature of necklace loops in our model is in many ways the opposite 
of what Matsuda predicted. We find that PBH formation is favoured at early times with potential 
dark matter relics forming later. This is indeed an unexpected result [20], though one which 
appears to follow naturally from the consideration of monopole/bead formation as a dynamical 
process rather than as the result of a separate phase transition prior to string formation. We 
argue that this is the correct approach to take when considering monopoles which form from 
extra-dimensional windings and a comparison of our results with field-theoretic string-monopole 
networks is given in section 5.6. 

^We assume here that the minimum radius of the loop is determined by the effective width of the string. 

^ds{ts) ~ ^M^B Unj) where ts is the formation time of the string network, tm is the time of monopole formation 
such that tiii < ts- 
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In addition to a renewed interest in the role of F / D-str'mg networks in cosmology there 
are still many in the physics community devoted to the study of field theoretic cosmic strings 
[211 [22l [23] . Whilst a stringy origin of the CMB perturbations has been ruled out by observation, 
the best A CDM model fit to the data suggests that these strings may contribute at the level 
of ~ 10%, 121] (for a review see also [25] ) making them extremely important objects to study. 
Unfortunately a best fit for (p, q) strings, necklaces or cycloops has not yet been investigated, 
at least at the field theory level |26j . 

However recent discoveries of dualities between gauge field strings and F/D-strings have 
raised the possibility of unifying these two approaches. Indeed, if string theory really is a theory 
of everything (TOE), and if we accept Quantum Field Theory (QFT) to be a valid low-energy 
approximation, we may also hope to find string theory analogues of all field theory phenomena. 
If therefore we expect topological defects, including strings, to arise generically in symmetry 
breaking processes, we must investigate the relationship between fundamental strings and field 
theory strings in much more detail. 

The aim of this paper is to consider a simple model of string necklaces in a well understood 
supergravity background using type IIB string theory, thereby extending the initial phenomeno- 
logical approach begun in [20] . Following the considerations above therefore we also consider the 
possible relation of these objects to field theoretic strings, restricting ourselves to generic consid- 
erations. Since the background yields an explicit form for the uplifting potential in terms of the 
number of extra-dimensional windings we may compute the bead mass precisely if this number 
is known for a string loop formed at any epoch. Following [TJ] we assume that the motion of a 
string in the compact space, prior to the chopping off of a loop, is random. This therefore allows 
us to estimate both the time-dependent bead mass and the average inter-bead distance. We 
find that the results depend the definition of the parameter uji, which gives the fraction of the 
total string length contained in the extra-dimensional windings. Two definitions are suggested: 
Identifying the inter-bead distance in the string picture with the correlation length in the field 
theory picture, we see that the first definition leads naturally to a scaling solution - similar to 
that for field theoretic strings but with a correlation length L(t) ~ much smaller than the 
horizon. The second leads to sub-scaling solution L{t) ~ t*. 

More importantly we are able to compute the PBH mass spectrum produced by collapsing 
necklaces. Since the mass of an individual necklace depends upon the structure of the internal 
manifold as well as the string tension, the resulting PBH spectrum yields information about the 
size of the of the extra dimensions and the warp factor. This in turn influences the background 
cosmic ray flux via the Hawking radiation of PBH's expiring at the present epoch. We are thus 
able to provide observational bounds on string theory parameters using measurements of the 
extra-galactic gamma-ray flux at lOOMeV from the EGRET experiment [27] . 

An interesting result is that both definitions of ui give the same qualitative behaviour with 
PBH formation occurring only over a limited time period in the early universe. However it must 
also be noted that the upper and lower limits of this window and hence the resulting bounds on 
the model parameters do vary significantly in each case. 

The layout of the paper is then as follows: In Section 2 we introduce the Klebanov-Strassler 
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geometry which is the supergravity background for the string embedding. In Section 3 we 
construct the world-volume action for string loops before analysing their stability in Section 4. 
Section 5 deals with the formation of these loops and their cosmological impact, focusing on the 
predictions for PBH abundance and Section 6 contains a brief discussion of our main results and 
suggestions for future work. 

We conclude this introduction with a note regarding terminology. The "necklaces" which 
are the subject of the present paper and of Matsuda's original work should not be confused with 
"necklaces" formed via other string-monopole interactions. For example, Leblond and Wyman 
|28j have shown that L>0-branes (monopoles) may be formed at junctions between strings in a 
(p, g)-string network. The resulting "necklaces" are of no relation to the ones considered here. 
Related work can be found in p9| [301 [3T| [32] . 



2 Klebanov-Strassler geometry 

The background we wish to consider is that of the Klebanov-Strassler (KS) throat |18[ 119] since 
it is one of the better understood backgrounds of the type IIB theory. Recall that conical singu- 
larities are the most generic kind of singularities arising within compactifications of IIB string 
theory on manifolds of SU{3) x SU{3) structure j33| . Since explicitly realistic compactifica- 
tions are difficult to construct, we will take a more phenomenological approach by considering 
(non-compact) conical backgrounds such as the conifold T^'^, that can be glued to a (conformal) 
Calabi-Yau manifold. Provided we work in a region far from this gluing, we can (locally) work 
with the conifold geometry without worrying too much about the precise details of the com- 
pactification mechanism. With this in mind we can regularise the singular conifold by allowing 
the to shrink to zero size. This is just the deformation of the conifold, which is topologically 
equivalent to the cotangent bundle over the three-sphere T*S^ where the has some minimal 
size. 

The solution of Klebanov-Strassler requires the introduction of both D3-branes and fractional 
Z)3-branes in this background [18], where the fluxes back-react on the geometry to create a 
warped throat. The RR-three form flux is threaded through the finite size three-sphere, whilst 
the NS-NS fiux wraps the dual (B) cycle as follows, 

[ f(3) = M, / ^ (2.1) 

A-K^a' J S3 47r^a' Jb 

The fiuxes generate a non-trivial warp factor for the full ten-dimensional solution which effec- 
tively measures the size of the blow-up contribution along the 5^ via 

h{T) = 2^/^{gsMa'fe-^/^I{T) 

I{t) = r ~ ^ (sinh(2x) - 2x) dx. (2.2) 

Jr sinh [x) 

For the solution we are considering, we must send r — t- 0, since r controls the radius of the 5^, and 
therefore the warping approaches a constant value ao = LimT-_5>oh(r), and the ten-dimensional 
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metric is well approximated by the following; 

ds^ = ali^^.ydx^dx" + R^{dilj^ + sin^ + sin^ Odcf)'^)) (2.3) 

with ~ C(°) ~ at the tip, and 

C(2) = —MaFN))uj2 (2.4) 
9s 

where uj2 is the volume form along the two-cycle parameterised by 9, (j), and F{^) ~ ^ — 
sin ^ COS'?/'. In this definition we are explicitly assuming the background gauge choice [S] 

6 = 9i = —02 

<P = (t>i = -(^2 (2.5) 

and is now the azimuthal angle on the with il^ € (0, 2-k\ as its fundamental domain. Note 
also that the are dimensionful coordinates whilst we define the radius of the three-sphere via, 

= bMgsa (2.6) 

where 6 is a numerical factor and gs is the string coupling, keeping the angular variables dimen- 
sionless. It is important to note that the dilaton, and therefore the string coupling, is constant 
in this background. 



In what follows we will assume that (|2.3p is representative of a large class of warped ge- 
ometries, without necessarily having explicit realisation in a fully UV complete construction. In 
particular this means that we should treat the warping parameter ao as a constant satisfying 
< ao < 1. For the (non-compact) KS solution we find that ao is related to the deformation 
parameter of the conifold e via the expression 

;4/3 

«o - ^ (2.7) 

where e is re-scaled to have the correct dimensions and is related to the deformation parameter 
e of the warped conifold. The basic point is that the warping varies like l/ggM where we 
are assuming the supergravity (SUGRA) limit gsM >> 1 so that, for fixed R, constraints on 
the warping can be interpreted as directly constraining the background geometry. We should 
also point out that there is nothing special about the deformed conifold solution. One could 
equally well use the resolved conifold where the is blown up instead \3A\ I35| . Indeed the 
resulting cosmic-string tension scales directly with the resolution parameter and is therefore 
highly constrained by observations. The resulting solution is then very similar to the model 
considered in [37] . 



Our choice for the background is also inspired in part by the AdS/CFT duality, since the 
KS geometry is known to be dual to an = 1 confining gauge theory. When viewed from 
this perspective, the strings are effectively the confining strings of the gauge theory with non- 
trivial wrapping. Whilst this is interesting in its own right, our motivation will be to model 
cosmic-necklaces rather than gauge theory necklaces. 
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3 String loops with non-trivial windings in internal space 



The action for both fundamental strings (F-strings) and L)l-branes (D-strings) - or in more 
general {p,q) notation (1,0) and (0,1) strings respectively - in the warped deformed conifold is 
simply the Nambu-Goto action with additional world-sheet flux, plus a possible Chern-Simons 
term (for the D-string), 

S = -Ti j dadt^/^ + Scs (3.1) 
where Ti is the string tensior|3 and X = detXab, 

Xab = lab + XFab (3.2) 

where a,b £ {0, a}, jab is the usual induced metric on the world-sheet, 

lab = GMN{X{a, t))daX^\a, t)dbX''{a, t) (3.3) 

and A = 27rZ^ . The simple form of the Lagrangian density \/—X arises because we are neglecting 
the coupling of the world-sheet to the NS-NS two-form field, which is vanishing in our back- 
ground (at least in the limit we are considering). The flux tensor is anti-symmetric such that 
A-Fqo = AF(j(j = and AFqo- = — AF^-o, AFoo- > 0. Note that we are absorbing the definition of 
the string coupling into the field strength tensor, since this allows us to identify the String frame 
with the Einstein frame. The Chern-Simons coupling is given by the integral of the pull-back of 
the RR two- form over the world-sheet, for the case of D-strings, 

Scs = Ti J dcrdtM a' (ip- sin ip cos i;)de A d(j) (3.4) 

However, for the sake of simplicity we choose to ignore this correction in the following analysis 
and instead concentrate solely on the Nambu-Goto component of the action. 



We wish to take the following general ansatz for the string embedding, 

= {t,r{t)sma,r{t)coscr,zo,T ^ 0,0,0,ip{(7,t),e{a,t),<j){a,t)) (3.5) 

where O<'i/'^27r,O<0<7r and < (/) < tt. This describes a circular string in the Minkowski 
directions which is wrapped smoothly over the in the internal space. Note that in this model 
the strings sit at the tip of the warped throat. We also make this choice for the sake of simplicity. 
However there is no a priori reason why we should make this assumption. The more general 
case would follow the lines of [SB]. Using the metric (|2.3p and the embedding ansatz (|3.5p the 
action then becomes 

S = -Tial J dadt^il - r'^)ir^ + a^'^R^s'^) - Oo^r^R^s^ - %^X^Fl - a^^R^s^s'^ + ao^^^(ss')^ 

(3.6) 

^Note that in the warped throat geometry we do not expect the tensions of the F and D-strings to be exactly 
equal, T(i,o) 7^ 2^(o,i)- The formula for the tension of a general (p,q) string, in the large p/large q limit, in the 

warped throat background, is T^p^q) ~ 2^ \j (^) + ■sin^ (^7)1 t)^^ it remains an open problem to calculate the 
tension for small p/small q. We may however assume that the order of magnitude of both the F and D-string 
tensions are set by the fundamental string scale, i.e. T^i,o) ~ "^(0,1) ~ a'^^ . Here we will use the term Ti to refer 
to either tension. 
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where we have introduced the sUghtly abusive notation, 

s2 = ^2 _^ sin^ '(/'(^^ + sin^ 6^"^) 

s'2 = ip'^ + sin^ il^ie'^ + sin^ e(t>'^) (3.7) 
{ss') = iPi;' + sin^ ^p{ee' + sin^ e(j)(p') 

where a dot or dash indicates differentiation with respect to t or cj respectively and we have 
chosen the gauge so as to identify the world-sheet time coordinate with the proper time in the 
Lorentz frame of the string loop, t. We have also included a non-zero 'electric' gauge field 
contribution for generality. For the D-string case this corresponds to the dissolving of F-string 
charge on the world-sheet and the strings are essentially superconducting. 



There are two constants of motion for this configuration, the total energy of the string H 
and the angular momentum I, due to the motion of the string in the internal dimensional^. We 
parameterise these conserved charges as follows [37] 

where L is the Lagrangian and are the canonical coordinates which are themselves functions 
of the world-volume coordinates through = {r{t),a,il^{a,t),9{a,t),(f){a,t)). Using ()3.6p and 
neglecting the Chern-Simons term, we find the following expressions for these conserved charges 

H = T, , al(,-^atR^s-) 

^(1 - r2)(r2 + 2i?2s/2) _ 2^2^2s2 _ a^^^X^Fl^ - a'^^ms'^s'^ + a^^m{ks'Y 



da 



^(1 - r2)(r2 + 2i?2g/2) _ a-2^2^2s2 _ a^^X^Fl - a^^R^s^s'^ + a^^R^{.ss')^ 



However it is more useful to re-write the Hamiltonian in canonical form using the momenta; 

dr J 













\dA 












\de. 



dL \ ^ ^9 / dL \ ^ , n / dL ^ ^ 



Ll = [^] , , ^^=(^) (3.9) 



This will simplify the form of our solutions for stable string windings with non-zero world-sheet 
fiux A-Fqo- 7^ 0. This is especially true in the static (/ = 0) case which we will consider in Section 

H 



After a long but straight forward calculation we find that the canonical form of the Hamil- 
tonian reduces to, 

®Note that although the string "rotates" around the S"^, no centripetal force is acting upon it. The internal 
(compact) dimensions are parameterised in terms of the angular variables 9 and 4> and so the motion through 
the 5*^ is measured in rad x \t\~^ , the units of angular velocity. 

®We have not here considered the more general case of a loop with non-trivial windings in the internal manifold 
which also rotates in Minkowski space. Such a scenario would require the action of a genuine centripetal force 
provided by the string effective tension. 



8 



where we have written 



TfalR^r'^ip + sin^ il;[6'{ipd' - dip') + sin^ 9ct)'{i^(j)' - ^ip')] 



f 2 



TfalR'^ sin^ sin^ 9r'^(l) + [^^'(0^/;' - ijj^)') + sin^ VlfA^" - 0<P')] 



In the next section we wih specify the ansatz ()3.5p completely by identifying the functions ip{(J, t), 
9{a,t) and 4){a,t). We will then use equation (j3.10p to determine the stability conditions in the 
static case by minimising the total energy of the string. Finally we see how this leads naturally 
to a concrete model of necklace loops in the warped throat scenario when these solutions are 
perturbed. 



4 Stability analysis for string loops in the static (/ = 0) case and 
necklace formation 



In the static case we assume that there is no motion in the compact dimensions and that the 
loop is neither expanding nor contracting in Minkowski space. Setting tp = 9 = (j) = Q and r = 
(or equivalently = = = and = 0) gives the static potential V (below) and I = 
as expected. 

V = j daal^Tf + ^^r^ + afR^s'^ (4.1) 

Note that even static strings of this kind may lose mass-energy over time due to the emission of 
gravitational radiation, causing them to shrink. However we refer here to the stability of extra 
dimensional-windings in the static gauge over an epoch in which the loop size r and energy 
V are roughly constant. The shrinking of necklace loops over cosmological time scales and its 
implications are dealt with explicitly in Section [5l 

We must now complete the ansatz (13.51] . Choosing il^{a,t), 9{a,t), and (j){(J,t) to be such 
that the angular winding is linear in a we see 

9{a,t) = nga + 9{t), (4.2) 
(pia,t) = n^a + 4>{t) 

^"This configuration corresponds to a situation in which the end point of the string is equally likely to move 
in each of the angular directions prior to the moment of loop formation. As we shall see in Section [S] when we 
consider the random walk regime, this assumption is well motivated from a physical point of view. We also note 
for future reference that in canonical coordinates, windings of this form do not wrap geodesies in the S*^. This 
has important implications for the tj-dependence of both H and I. 
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where n^, ng, ricf, ^ N and < cj < 27r then gives 



V = j daalT,^l + ^^r2 + 00^^2(4^2 + ^^^2 ^(^2 + ,^^2 ^^2)^ ^4 3) 

where and ^ are themselves functions of a and t according to ()4.2p . 

We see immediately that 1/ = ^(-0,0) indicating that the i;^)-direction is flat, whereas the 
■0 and 6 directions are "lifted" by the presence of a potential energy density V which is the 
integrand in ()4.3p . 

It is now possible to make a connection between wound strings at the tip of the conifold 
throat, and the cosmic necklaces predicted generically by Matsuda [13]. First we must deter- 
mine the conditions under which the wrappings described above are stable, which is done by 
minimising the total energy V. Perturbations of this configuration are then seen to give rise to 
beads, whose mass can be estimated from the functional form of the potential. 

However rather than simply minimising the potential, it is useful (and easier) at this point 
to develop a physical intuition for the string configuration. Treating -0 and 9 as independent 
variables, we may sketch the lifting potential (strictly speaking the lifting potential energy 
density, but from here on these two terms will be used interchangeably) which the string "sees" 
throughout the whole S^. This is done by plotting the integrand V in ()4.3p . The result for the 
principle range (0 < '0 ^ 27r, < < vr) is shown in Figure 1 for the purposes of illustration. 
The critical points and associated field masses may be found by diagonalising the corresponding 




Figure 1: Plot of the static potential V setting the tension pre- factor term to unity and also 
= 1 and a^^^R"^ = 1. 



Hessian matrix. We verify that V has equal local maxima at positions {ip = 2 ' ^ ~ 
(2m+i)7r saddle points of equal magnitude at {^p = (2"+^)^ ^ q = rnn) and flat directions which 
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are also local minima (i.e. "troughs" not "ridges" or points of inflection) given by {ip = mT,6), 
where m,n £ Z. This gives rise to two local maxima, two saddle points and two flat directions 
in the principle range. The associated field masses are given below. 

m tp = (^"+1)^ ^ g — (2m+i)7r ^jQpg^^ maxima) 



mi = ,m'^ = , (4.4) 

r2 + ao 'i?2(4n2 + + n^) ^r^ + ^R^i4nl +nj + n^) 

^ _ (2n+i)7r ^ Q _ (saddle points) 

r2 + ao ^ii2(4^2 ^ ^2) ^^2 + ^^2(4^2 ^ ,^2) 

V' = nn, 9 (flat directions) 

me = 0, m^ = — , ^ (4.6) 

'■2 + 4ao'i?2n2 

It is now intuitively clear that minimal energy configurations correspond to strings wrapping flat 
directions in the two-dimensional sub-manifold described by ip and 9. Physically this corresponds 
to strings wrapping some point in the 5^ which is uniquely determined by the condition tp = nvr 
for some n € Z. This may be seen from the metric (|2.3|) and is the reason why ng and do not 
contribute to the total energy. Even though, technically, ng,n^ > windings around points have 
zero length and are not physically meaningful. This result is precisely what we should expect, 
since the minimum energy configuration ought to correspond to a situation in which the string 
has effectively zero length involved in extra-dimensional windings. 



Indeed, we may verify this by substituting tp = nvr and setting n2 = in (|4.3p that the total 
energy of a string wrapping flat directions in the potential V is given by 

V = 2-KalrTi (4.7) 

which is simply the rest mass of a string loop with radius r in warped Minkowski space. 



This is certainly a very complicated way of verifying an intuitively obvious result but the 
above analysis will prove useful when we consider perturbations which result in different sections 
of a single string lying along equivalent local minima. That is, when the string interpolates 
between degenerate minima (flat directions) in the {^,9) sub-manifold resulting in the formation 
of beads. 



Note that the flat directions in Figure 1 effectively deflne degenerate vacuum states for the 
string. For this reason we propose identifying the inter-bead distance in the string picture with 
the correlation length of field-theoretic strings. This forms the basis for identifying the field- 
theoretic parameter 7 (which defines the correlation length as a fraction of the horizon distance) 
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with the parameters which define the KS geometry in the following section. Note also that, 
although we refer to the bead-forming states as perturbations from the minimum energy (zero 
winding) configuration, this is true only in an energetic sense. Such configurations differ locally 
from the minimum energy (zero winding) configuration only in the vicinity of a bead. 

Indeed it is unclear whether it is possible to create a necklace from a standard F or D-string 
loop via physical perturbations of a section of string after the formation of the loop. This is 
the same as asking whether it is meaningful for an open string section (which may itself form 
part of a string loop or part of a string section connected to the network) to contain fractional 
windings which result in the formation of beads. 

It is subtle point but in his generic argument Matsuda [13] predicted that, due to the presence 
of the lifting potential, integer windings would not be necessary for stability. In principle 
this should allow for the formation of beads in open string sections due to local fractional 
windings. Though there is no reason why this could not happen generically, the possibility is 
not explicitly realised in our model. As discussed in the following section, we assume that any 
small-scale structure due to local partial windings will quickly disappear due to the annihilation 
of beads/anti-beads. 

In our model the existence of the potential barrier which stabilises the windings after the 
loop chops off from the network (and which allows for the formation of beads) depends upon 
the presence of non-zero integer windings at the moment of loop formation. What is more, this 
ought to be true generically in this class of models since the form of the lifting potential in 
the extra dimensions ought always to be determined by the parameters that characterise the 
internal windings. 

Physically what happens is the following; before the loop chops off from the network the 
string is free to move and create windings in the compact dimensions. If the internal manifold 
were not simply connected these would become topologically trapped resulting in the formation 
of cycloops jl4j . However as the 5^ is simply connected, these windings cannot be topologically 
stabilised. Instead they are stabilised by the presence of the lifting potential V, which is itself a 
function of the number of windings. 

We may imagine that, at the moment of loop formation the string ansatz is described accu- 
rately by ()4.2p which is sketched in Figure 2a. However as soon as the string chops off from the 
network to form a loop, the total energy (V) is no longer minimised for such a smoothly varying 
configuration. Although windings may continue to vary smoothly in the (^-direction they will 
then adopt a step-like configuration in the (V', 0) sub-manifold which is illustrated in Figure 2b. 

Note that this is an approximation. Technically if the string configuration (i.e. the ansatz) 
changes, then the form of the lifting potential also changes. This results in a complicated 
iterative process with complicated string dynamics before the eventual formation of a steady 
state. 

However if we consider V to be approximately constant, we see that the integral V = J daV 
is minimised precisely for the step-like configuration shown in Figure 2b in which the string 
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interpolates between degenerate minima by crossing the potential barrier perpendicularly at its 
lowest point i.e. at a saddle point ij: = IHlLtllZL^ Q = f^Tr. This in turn may be thought of 
energetically as a series of small perturbations away from the genuine minimum energy con- 
figuration of zero-windings (described above) which is equivalent to a string wrapping a single 
flat direction in the {'ip,0) sub-manifold. Looking at it from this perspective helps us to justify 
our original assumption that V remains approximately constant, so long as we remember to set 
= in ()4.3p . This leads to an apparent contradiction. We may assume, at the moment of 




Kig2a rii2b 



Figure 2: Figure 2a shows strings winding smoothly across the effective potential V on the (9, <j)) 
sub-manifold. Figure 2b shows step-like configurations which minimises the total string energy 
within V. 

loop formation, that n ~ ~ ~ > 0. After this point windings in the (/i-direction will 
still be able to vary smoothly, though they will not be stable and may in principle contract until 
they are point-like. This possibility is discussed in detail in section 5.6. 

However, so long as > 0, windings in the 6) sub-manifold will relax into the step-like 
configuration shown in Figure 2b resulting in the formation of approximately 2n^ beads. This 
forces us to view the configuration as a series of small perturbations (in fact a series of 2n^ 
perturbations) away from a configuration described hy tp = mr and = 0. We must therefore 
continue to regard > (the number of windings in the ^/^-direction at the moment of loop 
formation) as physically meaningful when estimating the number of beads, but we must regard 
it as approximately zero when estimating the bead mass. This is a little strange but it is quite 
consistent with the physical picture we have been sketching. 

To obtain our estimate for the bead mass (Mf,) therefore we first Taylor expand ()4.3p with 
set equal to zero. Assuming agr >> nR at the moment of loop formation, which corresponds 
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to a large loop in the Minkowski directions, for all times we obtain the approximate expression, 



We then move from the smooth windings picture shown in Figure 2a to the step-like configuration 
in Figure 2b by setting tjj = mn globally, except in the vicinity of a bead where = mvr, dij: = da 
and mTT <il}<{m + 1)-k. Our expression for the total energy V is then, 

I / 1 i?2 r-<pj={m+l)iT \ 

^ ~ W 1 + ( 2'nalrTr + 2n^ x -Ti — W il^d^ \ 



~ Wl + ^ [2^alrT^ + 2n^ x jT^—^'j . (4.9) 

As the first term is simply the rest mass of the string in the Minkowski directions, the second 
term corresponds to the rest mass of = 2n^ beads. Finally, setting ~ ~ n explicitly, 
we have the following estimates for both the initial number of beads N and the bead mass Mf,, 

iV ~ 2n 




(4.10) 

Note that the bead mass is inversely proportional to r. This makes sense on the assumption 
that the total length of the string remains constant. So if there is less string length involved 
in internal windings (resulting in less massive beads) more length is added to the ordinary 
four dimensional part of the loop. The converse is also true. Thus the first expression for the 
mass increases in magnitude as the second decreases and vice-versa. Note also that the M;, is 
effectively quantised in terms of the number of windings present in the ^-direction. 



Let us now also consider the full expression for the bead mass. One can integrate the 
potential over the above range and the result is roughly the sum of the bead mass and the rest 
mass associated with an unwound string (with n = 0). We can then re- write the bead mass 
including all the higher order terms, with an appropriate normalisation to yield 

M. ^ 2agT,.^l + ^ (ElhpticE (^) - ^) (4.11) 

which is valid for non-zero winding number n. One can easily check that in the limit where 
nR « a^r the mass reduces exactly to the one in (|4.10|) . 



5 Cosmological implications of necklace loops 



We now investigate the cosmological implications of necklace loops based on the assumption that 
they retain their necklace structure after formation. As we have already seen, the ^-dependence 



14 



of the lifting potential depends on the presence of windings in the (/)-direction (n^ > 0) and its 
■0-dependence requires the presence of windings in the ^-direction {ng > 0). However since the 
(^-direction is flat, windings along this direction are unstable. If these windings contract this 
effectively flattens the ^-direction leaving the ^-windings free to contract as well. This in turn 
flattens the ■0-direction and necklace structure disappears as the bead mass (which comes from 
extra-dimensional windings) is converted into the ordinary rest mass of a loop in Minkowski 
space. We therefore see that, due to the presence of single flat direction in the 5^, the entire 
necklace structure of the loop is unstable and may unravel in time. 



We will consider this possibility later in the present section where we will introduce a time- 
dependent model for the number of windings. For the moment we will we will consider n ~ 
n^j} ^ ne ^ to be roughly constant. We can now use Matsuda's original assumption that the 
four-dimensional part of the loop loses mass-energy via the emission of gravitational radiation in 
Minkowski space (just like an ordinary cosmic string in four dimensions) but that the bead-mass 
which is formed from the winding of the string in internal space is unaffected by this process. 



To investigate the cosmological implications of necklace loops we must therefore modify 
equation (|4.9|) by inserting a time-dependent radius r(t,tj) into the first term of the expansion 
(the loop mass in Minkowski space) and the initial radius r{ti) into the second term (the bead 
mass). The time dependent radius of a shrinking loop in warped Minkowski space is given by, 

aor{t,ti) = ao{ati - TgTi{t - U)) (5.1) 

where ti is the time of loop formation. Note that t > ti is the cosmic time coordinate, F is a 
measure of the rate of energy loss due to the emission of gravitational radiation, Q is Newton's 
constant (which is determined by the volume of the S^) and < a^a < determines the 
characteristic initial loop radius as a fraction of the horizon, dn ~ aoct The initial loop radius 
as a function of ti is then 

aor{ti) = aoati. (5.2) 

The bead mass now explicitly depends on the time of loop formation, ti via ()5.2p . However we 
also expect the initial number of beads, N, will also depend in some way on ti. We therefore 
need some way of estimating the initial number of windings in each direction n = n{ti). 

Following Shellard and Avgoustidis |12j(and taking into account the warp factor oq) we use 
the random walk regime to estimate the initial number of windings present in each angular 
direction, n(ij), so tha10 

^aoaujieiti 

n{ti) = (5.3) 



Although in general a random walk will not give rise microscopically to completely smooth windings, as 
described by the embedding ansatz, it should on average produce something similar from a macroscopic point 
of view. Any extra beads formed by microscopic structure would quickly annihilate one another if we assume 
that they are free to move around the shrinking loop in a random walk. This random motion occurs when 
small sections of the string momentarily acquire enough energy to jump between adjacent minima in the effective 
potential, causing the bead to move from a four-dimensional perspective. In his original paper [IS] Matsuda 
predicted that, based on the idea of a random walk, approximately \/n{ti) beads/anti-beads created in pair 
formation would survive until late times. However we contest that this is valid only for a static loop. For a 
shrinking loop it seems clear that all bead-anti-bead pairs will eventually collide and annihilate, at least if the 
minimum radius of the string (~ string thickness) is comparable to the initial spacing. This implies that only 
beads formed from net windings will contribute to the mass of any PBH's/DM relics eventually created. Further 
discussion of this point is given in Section 6. 
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where < < 1 is the fraction of the total string length 1 (not to be confused with the angular 
momentum) contained in the windings and €i is the step length. This definition tells us that 
string networks will tend to form once the correlation length becomes larger than the scale of 
the internal dimensions. 



The parameter uJi may be defined as the ratio of the string length in the internal dimensions 
to the total length of the string in all dimensions. Following |12j . but using the notation defined 
earlier, the concise definition of uJi is given by 




^i = \l I {-2-^;^ ]^da/ I Vda (5.4) 



where V(cr) denotes the integrand in equation ()4.3p . In order to simplify the above expression 
we then make the following approximations, 




/ n'^R'^ /(4 + sin^(no-) + sin^(nc7))dc7 

where in the second step we have ignored the numerical coefficient in front of the term, and 
again used assumption that n ~ ~ rie ~ n,^. Finally we note that it is straightforward to 
show that 

y"(4 + W{na) + sm'^{na))da ~ 0(10) (5.6) 

for all possible values of n[ti). Therefore we see that loi can be well approximated by the following 
function 

nR 

up to various numerical factors which would have appeared in the above expression if the integral 
in (j5.4p been performed in full. It is clear however, that w/ has the correct functional dependence 
on both n and r. This means that, whatever time-dependent models we use for r(tj), n(tj) and 
ij^iiti) - they must satisfy ()5.7p for consistency. Note that at this stage we are not imposing any 
additional conditions on a^r and therefore loi is just a parameter of the theory. 



If we then use (|5.2p as our model for r(tj) and substitute oj; from (j5.7p into our previous 
expression for n{ti) (which also uses the approximation r(ti) ~ atj) we find a cubic equation in 
the variable (nR)'^ ^ 

{nRf + al(atif{nRf - al{atifef{nRf = (5.8) 
This has the trivial solution (nR)^ = (no windings) and the more physically relevant one, 

'-)-^(-/^) 

although the reality of n requires us to take the positive sign before the square root. The time 
dependence of n{ti) in the expression above is then consistent with the definition of uji in ()5.7p . 
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Next we move on to consider to size of the step length e;. The maximum velocity of the 
string in the compact dimensions is c (where c = 1 in natural units), which corresponds to a 
step length per unit time (At = Oq = Oq^Vo') of ei = a^^ls = aQ^Vo/. El However only 
the endpoints of the string at the horizon move at the speed of light. For two points on the 
string within the horizon, separated by a distance d ~ act, the relative velocity between them 
is I) ~ ac which corresponds to an effective step length of, 

e; ~ aaQ^Vo/ (5.10) 

Note that this definition also implies that n(tj) oc a as we would expect. This ensures that 
the number of beads on a long string, which stretches across the entire horizon riHoriti), is 
independent of a. In fact the assumption of a constant step length e; at all points along the 
string is problematic, as this leads to a measure of nHor{ti) which is proportional to 1/^/a. 



Using (j5.10p and ()5.9p the condition for bead formation, n{ti)^ > 1, is then equivalent to, 

ti > : (5.11) 

and reality of the above solution translates into the constraint 

2 / 

n a a 

«o < -J^- (5-12) 

Now strictly speaking, Qq is fixed by the ratio of the fluxes arising in a full string compactification 
and is therefore highly sensitive to the magnitude of string coupling, the string scale and the flux 
parameters. In the non-compact case which we are considering we also require that Oq < 1 in 
order to ensure that the solution is warped. Therefore the existence condition for beads appears 
to impose a strict upper bound on the value of Oq which in turn imposes a strong constraint 
on the deformation parameter of the conifold geometry. We also note that if Oq = this 
allows necklaces to form only at infinity, and so only the strict inequality in ()5.12p is physically 
meaningful. 



One can now substitute the solution for nR into the generalised mass function in (|4.1ip . If 
we define the following function 



4a' 

G{U) = -l + Jl + ^ (5.13) 
then we see that there is a remarkable cancellation of terms and we are left with 




Nh ~ 2alT,aU^l + ^ (^EUipticE ( i\l ^ ) - J ) (5.14) 



'^■^Here we have used the fact that the standard string tension Ti ~ is given by the fundamental string 

mass divided by the fundamental string length Ti = -y^ ~ ^'^^L- ~ a'~^ together with the fact that the effective 
tension of the string in the warped throat is Ti ~ OqTi. We then treat this as the division of the 'warped' 
string mass riis ~ aov^ by the 'warped' string scale le ~ a^^y/a' so that the velocity of light is given by 



V a' 



Is V 
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which is not a monotonic function of time. For vanishingly smah ti the mass is increasing until 
it reaches a maximum value, before decreasing monotonically. This leads to the interesting 
possibility that one can have sustained PBH formation during a small window where the mass 
function is greater than the critical value required to produce a gravitational radius larger than 
the string width. This possibility is dealt with explicitly in section 5.4. 



We also note that the explicit form of uji{ti) is given by, 



.,~^(-l + ,/l + i;^l=^G(t.) (5.15) 



which tends to unity as tj — )■ and zero as — )• oo as we would expect. 



To aid visualisation in the discussions that follow we plot, in Figure 3, the functions Mh(tj), 
n{ti) and uJi{ti) and the inter-bead distance d{ti) for fixed values of our model parameters to 
illustrate the qualitative behaviour of each. What is clear from the plots is that, assuming no 




Figure 3: Plots of the mass of an individual bead as a function of time ti, the number of 
beads n{ti), uji{ti) and the average inter-bead distance d{ti). In all plots model parameters have 
been set such that qq = 0.1, a' = 1, R = 5, a = 0.8 and TiJl + = 1. 



gravitational emission, there are initially few beads with large mass. However as time increases, 
the number density of beads gets larger and their corresponding mass decreases. The rate of 
decrease is in fact much larger than the increase in the number density. Therefore one expects 
that eventually the system will become dominated by a large number of (almost) massless beads. 
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5.1 Late time formation - the scaling regime 



Using the fact that nR is fixed by the parameters of the theory, let us consider the asymptotics 
to better understand the physics at early and late epochs. Let us initially focus on the late time 
regime where 



ti » (5.16) 



which, using (|5.9p . implies that the initial number of beads is constant 

2 / 

n a a 

Note that the condition for bead formation n{ti)'^ > 1 simply reproduces ()5.12p . ensuring the 
reality of i^'"*"') . the minimum time at which necklaces formation begins - and is consistent with 
our expectations. We are then able to estimate the average inter-bead distance for necklaces 
formed at late times 

d{U) - ^ - (5-18) 

Identifying this with the correlation length along the string, L{ti), as mentioned in the previous 
section, we see that our solution corresponds to a scaling regime in which we may identify the 
scaling parameter 7 with the Klebanov-Strassler parameters, i.e. 

d{ti) ~ L{ti) ~ aojti (5.19) 

implying that, 

7 - ^ (5.20) 

Compare this with the result of field theory calculations, where the initial number of beads on 
a loop of radius r(ti) in the scaling regime is given by |16j . 

, , r(tj) a 
^[k) 7 

Hence we see that the condition for bead formation in the string picture (|5.12|) is equivalent 
to the condition for bead formation in the scaling regime of the field theory, 7 < ^a. Both 
conditions ensure that N > 2, implying that beads come in pairs. This is also consistent 
with the identification of the inter-bead distance with degenerate vacuum states, or correlation 
distance, along the string. Thus we see that a scaling solution arises naturally at late times 
in the string picture allowing us to identify various string parameters with their field theory 
counterparts. 

Using the condition on 7 and a which is consistent with bead formation i.e. < 7 < ^a, this 
simply recovers our previous constraint on gqU (|5.12p . However by setting a to the maximum 
value allowed by causality a = 1 we may place a maximum upper bound on the radius of the 

-1 



R < a^Wa'. (5.22) 
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This simply means that beads/ windings wih not form unless the radius of the is smaller than 
the warped string scale, irrespective of the value of a. It also implies 



«o < TTl— (5-23) 

in order for windings to form. Note also that setting oc a ensures that n{ti) oc a even if uji is 
left as a free parameter as in [12]. Otherwise we have that n(tj) cx y/a which implies that the 
number of windings per Hubble volume is proportional to \foL which seems unphysical. 

We now ask what happens to the time dependence of the mass function in the late time 
limit. The argument of the EllipticE term becomes small and therefore we may series expand 
the mass term to obtain 



n2 vra^a't / 3 



^^~T'i^/l + ;^--3- 1---^ + ...) (5.24) 



and therefore we see that the mass evolves inversely proportional to tj and that M;, — )■ as 
— >• 00. This implies that any necklace structure should disappear at sufficiently late times. 
In other words it is unlikely that any necklaces formed at late times would be distinguishable 
from ordinary string loops. One should quantify this by noting that the mass is also inversely 
proportional to the scale of the warping, therefore for highly warped throats the mass will be 
constant over a much larger time scale. 



5.1.1 Early time formation 



Returning now to the cubic solution for (ni?)^, let us consider early time formation subject to 
the condition that ti < . This gives us 



B? 2\R 

to second order, which also implies that, 

1 

2Va' 



ioiiU) 1 - ^^alu (5.26) 



from the definition of uji. Clearly therefore c<;; — t- 1 as — t- 0, that is as aQr{ti) — t- 0{nR). 
Utilising the bead-formation condition then allows us to place a bound on the bead formation 
time via 

, , , <U< (5.27) 

The average distance between beads, should necklaces begin to fornix. . may then be approxi- 
mated by 



^^Note that we have assumed t^J^'^'^^ ^ — < ^/p" which is equivalent to assuming aoR < O.Olav^ 
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which, unhke the the late time approximation, does not correspond to any known regime in the 
field theory if we continue to identify d{ti) with the correlation distance L{ti). However we note 
that the second term is sub-dominant at very early times when ti < suggesting that it may 

be reasonable to keep only first order terms in the expansion especially if << Such a 

very early time approximation may correspond to a damping regime in the field theory picture. 
The above equation would then represent an intermediate regime where damped and scaling 
solutions can be joined together. Assuming now that ti « and keeping only the first 

order term in ()5.25p . we see that the condition for bead formation r?{ti) > 1 is then equivalent 
to, 

U > (5.29) 

This is consistent with our previous estimates and shows that bead formation will occur in the 
very early time regime0only when a^B? « a^a' . With this in mind we see that the inter-bead 
distance is set by the leading factor in ()5.28p which is consistent with the field theory picture 
at early times during the damping regime, where we expect the correlation distance L{ti) to be 
given by a power law solution of the form [16] , 



L{U)^tlt] (5.30) 
with t(i corresponding to the characteristic damping time of small scale oscillations on the string. 



However a damping regime is usually obtained by considering collisions of the string with an 
external plasma. In this case the damping term comes from the internal dynamics of the model, 
suggesting that the inertia of the beads (when M^, is large) is sufficient to cause the correlation 
length to scale as L(ti) ~ y/U at very early times. 



Furthermore unless the warping is extraordinarily large, the time-scales over which this effect 
takes place are likely to be insignificant in cosmological terms. As a future amendment to the 
current work therefore it will be useful to impose an external damping regime and to study 
the effect of the interaction of the windings with the external plasma. Naively we may expect 
collisions of the string with particles in the compact space to inhibit the formation of windings 
resulting in the delayed on-set of a scaling regime. This to is likely to mirror the field theory 
case, though further investigation is needed to establish whether the correlation distance scales 
according to (|5.30p . 



Once again, for the sake of completeness, we consider how the mass function changes as a 
function of time in this epoch. The elliptic integral is actually divergent in this limit, however 
we can consider the leading order divergence which will dominate the spectrum. The resulting 
expression for the mass function will then behave like 



/ tt2 

Mb ~ 2Ti vW^W 1 + ^ {o^Uf'^ (5.31) 



^'^Note that the condition^ ~ „^>— < imphes that a^R < \ploL^/a! which is automatically satisfied 

by ^M- 
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which one can see will tend to zero as ti goes to zero. This is within the regime where we expect 
PBH formation to occur, since the mass of the necklace plus beads is steadily increasing (as a 
function of time) during this regime. 



5.2 PBH formation 



We now calculate the contribution to the PBH mass spectrum from collapsing necklace loops, 
based on the assumption that loops which chop off from the string network retain their necklace 
structure indefinitely. We will initially take the number of beads per loop to be constant from 
the time of formation. 

The minimum radius to which a contracting loop may shrink, 6i, is limited by the string 
width - which we assume to be comparable to the inverse of the symmetry breaking scale r]s 
I20j. i.e. 6i ~ rjj^ El. The condition for gravitational collapse then becomes 

Rs > 1]-^ (5.32) 

where Rs is the Schwarzschild radius of the loop. We may estimate the Schwarzschild radius of 
a necklace loop using the spherically symmetric approximation 

Rs{k) ~ 2QMT{ti) (5.33) 

where Q is the modified Newton's constant and Mxiti) is the total mass of the necklace. A 
necklace formed at time ti will therefore collapse to form a PBH if, 

Mriti) > ^ (5.34) 

Assuming that 6i is small, the bead mass will provide the dominant contribution to Mxiti) so 
that 




(aUfy^) [^ElhpticE l^i^J^j - (5.35) 

which is not a monotonic function. We can approximate the solution at early and late time 
regimes respectively using the following 



n2 vr a^a't 1 / 3a' 



^ + r2A22 alR U V 16a^t2 ) (^-^S) 



It can be shown that, for Abehan field theories, the tension of a string is proportional to the square of the 
symmetry breaking scale, ~ if . Our Abelian string with the same assumption gives 77s ~ \/q' However we 
will leave t)s as a free parameter in the discussion that follows. 
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Here we note that Mt is linearly increasing with a^ti at very early times, whilst in the late 
time epoch we find that the scaling is with l/{a^ti). Thus the sensitivity to the warp factor is 
most pronounced at late times, since a vanishingly small value of oq means that the total mass 
remains larger over a wider time-scale. Ultimately however, the total mass will tend to zero 
asymptotically. The complete mass function Mxiti), using the definition of the elliptic integral, 
is sketched in Figure 4 to illustrate the time dependence. In Figure 5 this is shown together with 
the early and late time time approximations. Without precise values for the model parameters 

1.0 n — ^ — ^ — ^ — I — ^ — ^ — ^ — I — ^ — ^ — ^ — I — ^ — ^ — ^ — I — ^ — ^ — ^ — n 
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Figure 4: The total mass Mt contained within the beads of a necklace as a function of ti. The 
model parameters have been fixed so that oq = 0.1, a = 0.8, a' = 1, Ti ^1 + ^^^^ = 1 and 
R = 5. 



it is difficult to determine the accuracy (as a percentage estimate) of either the early or late 
time approximations. However it is interesting to note that the late time approximation also 
qualitatively captures the behaviour of Mj'{ti) in the early time regime. In fact the peaks of the 
two functions are in approximately the same position, although the peak of the approximation 
is roughly twice as high as the peak of the true curve. 



However it is also clear that the peak of the full mass function, which is potentially the region 
of greatest interest with respect to the formation of PBH's, lies between the regions in which 
either the early or late time approximations are valid. As we shall see, in order to calculate the 
PBH spectrum analytically in terms of our model parameters, we must integrate over the region 
of the curve (|5.35p which satisfies ()5.34p . This is not possible in general, as any such region (if it 
exists) will certainly include the peak itself. We must therefore find some way of approximating 
the elliptic integral in this crucial region. 



Whilst this can be done in many ways, we will use the following method. To begin with, 
we estimate the maximum height of the full mass function by utilising the fact that the peak 
of the late time approximation lies at roughly the same value of tj. Differentiating the second 
equation in (j5.36p and solving the resulting expression dHx/dti = gives, 

^{peak) 3 

4^^ 



tr'> ~ (5.37) 



^0 
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Figure 5: The total mass Mt contained within the beads of a necklace as a function of tj, together 
with the early and late time approximations given in equation ()5.36p . The same parameter 
choices were made as in Fig. 4. The blue (green) curves are the early (late) time approximations 
respectively. Note that the late time approximation seems to over-estimate the maximum mass. 
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Substituting this back into ()5.35p then gives, 



^(peafc) ^ ^Titt^'^) ~ 5.47 Ti^l + (5.38) 

We then note that substituting the fuh expression for n(tj) ()5.9p into ()4.10p and using Mxiti) ~ 
N{ti)Mf){ti) from ()5.35p also gives a function which shows the same quahtative behaviour as the 
true Mxiti) curve (as we would expect). Taking this approach corresponds to expanding the 
elliptic integral only to first order in n{ti), but keeping the full time dependence of this function 
(which is valid even at early times) as opposed to keeping higher order terms in the expansion 
of the elliptic integral and using the late time approximation n(tj) ~ a^fa' /a^R as in (j5.36p . 
The resulting expression for Mxiti) is, 




Mr(t,) ~ ^T,y 1 + ^-1 + y 1 + -^J (5.39) 

Although this is still highly inaccurate within the region of the peak, a function of this form 
may be used to capture the behaviour of the full expression right down to all but the earliest 
times (where the early time expansion above must again be used)0. We may then fit a curve 
of the form. 



where A and C are free parameters, to the true Mxiti) curve by demanding the approximate 
function satisfies; 




• to pass through the point (p'^'^*^) ^ m^^""*^^ ) = (4f|,5.47 Ti ^1 + ^^^^ slo'j ) 

• to have the same asymptotic behaviour as as the full elliptic integral. This is equivalent 
to the requirement that A = -^C~'2. 



These two condition are sufficient to fix A and C uniquely giving, 

A 0.368657 

C 5.25648 (5.41) 

The resulting fit is remarkably good both at late times and in the vicinity of the peak and is 
shown for a selection of parameter values in Figure 6 below. Only at very early times does the 
approximation appear to break down, which becomes clear if we 'zoom in' near tj = as shown 
in Figure 7. However, if we are required to integrate over a time range that includes a region in 
which MTapproxiii) given in ()5.40p becomes invalid, we may always split the resulting integral 
into two parts. It would then be necessary to integrate, using the early time approximation in 
()5.36p and expanded to arbitrary order, over some range < ti < tE (where Ie marks the time 
at which the approximation starts to break down) - whilst using the expression (|5.40p / (|5.4ip 
over the remainder of the range tE ^ ti < tp- 
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Figure 6: Plot of the exact Necklace mass Mt (red curve) and the approximate fit function 
Mxapprox (green curve)as a function of U for oq = 0.1, 0.05, 0.01, 0.001 (top left to bottom right). 

In all the plots a' = 1, R = 5, a = 0.8 and TiJl + = 1. 
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Starting with (j5.32p we may now estimate the range of time over which necklaces will form, 
that will eventually collapse to produce black holes i.e necklaces with sufficiently high mass 
contained in their beads to produce a Schwarzschild radius greater than the string width. From 
our previous discussion we see that the mass function peaks at earlier times, therefore favouring 
PBH formation in this regime - that is, from the collapse of small necklace loops with a relatively 
low number of high mass beads. 

In fact it is clear from N{ti), Mi,{ti) in Figure 3, and especially from the plot Mxiti) ~ 
N[ti)Mb{ti) in Figure 4, that the increase in bead number at late times is insufficient to com- 
pensate for the reduction in the mass of individual beads. Similarly at very early times, there 
are an insufficient number of extra-dimensional windings to produce enough beads to create a 
large Schwarzschild radius, even though the individual bead mass is high. This results in the 
formation of a window in the early universe, during which necklaces suitable for the formation 
of PBH's are produced. 

The limits of this window, let us call them t^^jj, may be estimated by solving equation 
(j5.34p with the appropriate expression for Mxiti). Let us now assume that the energy scale 
of the symmetry breaking process which gives rise to the (p, q) string network is equal to the 
un- warped string energy scale, 

rjs ~ Vc/^^ (5.42) 

'^^It is more useful in this respect to keep the full time dependence in n{ti) while expanding to only first order 
in n than it is to expand to two or more orders in the argument of the EUipticE for large ti. 
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and consider the limits which arise from inserting ()5.40p with the values of A and C given 
above. The resulting time bounds can be solved for analytically but yield rather complicated 
expressions which we do not repeat here. 



These analytic expressions are valid so long as the equation ()5.40p provides a relatively good 
approximation to the true Mr(tj) curve. However the validity of the lower bound, tp^^j, should 
always be checked (for a given set of parameter values) by plotting the full mass function together 
with the PBH formation bound, 

' (5.43) 



2Q 2 G(io) 

where we have used the fact that the effective four-dimensional gravitational coupling is related 
to the ten dimensional one through the relation 

where Xq is the internal manifold. We can approximate the coupling via 

(5.45) 



although we will typically assume that G'(io) ~ C'(l) (in units where a' = 1) for simplicity. In 
the event that this forms a poor estimate of the true tp^jj it is necessary to expand the early 
time approximation given in (j5.36p to as many orders as required to reach an accurate value 
and to re-solve ()5.34p . 

For the sake of completeness we note that the early time expansion given to first order, as 
in ()5.36p . remains reasonably accurate up to, 



tE-^^ (5.46) 

beyond which the function ()5.40p is undoubtedly valid. The resulting estimate of tp^j^ (where 
tpBH <iE^ ^) is, 

t'pBH / (5-47) 

In order to be safe therefore we may choose to split up any integral involving the the necklace 
mass MT{ti) into two pieces: the first using the early time approximation ()5.36p between the 
limits tj,pfj < ti < Ie and the second part using ()5.40p between tE < ti < t^^fj where Ie is 
given by (|5.46p . 

An estimate for tpp^j can also be obtained using the late time expansion for MxiU) and this 
yields the simple expression:- 
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The next step towards calculating the necklace contribution to the present day PBH spectrum 
is to estimate the time at which a collapsing loop reaches its minimum radius 6i ~ 7?7^ ~ Vo'- 
This depends on the time the loop formed ti via 

We expect, from our approximations, that 6i will be small and therefore the second term will be 
sub-dominant. This implies that the corresponding time range over which PBH's actually form 
from the collapse of these loops is given by t'p^fj < t'{ti) < t'^^fj where 

(l + ^)ipijH- (5-50) 

The initial mass of a black hole forming at t'[ti) is therefore equal to Mxiti), which can be 
compared with the mass of black holed formed from ordinary string loops - MpBH{ti) ~ jJ'Odi-, 
which yields the present day mass spectrum |38j 

^"^;fi^'^ oc M-^-\ (5.51) 
dM ^ ' 

Clearly the mass spectrum for black holes formed via the collapse of necklaces is far more 
complicated, although the relevant calculation of the present day spectrum proceeds in a similar 
fashion. Firstly we identify 

MpBH{t'{t.i)) = Mriti) (5.52) 

Assuming also that the PBH formation window lies in the radiation dominated epoch, the 
number density of string loops with initial length r(tj) = ati which chop off from the network 
at time ti is given by, 

n(r(ii),ti) ~ ^ — r = -fr (5-53) 

where Vr is the number of long strings per Hubble volume. Now in field theoretic models, 

Ur = g-f~^c (5.54) 

where ^ is a Lorentz factor, 7 gives the correlation length as a fraction of the horizon (which 
we previously identified with the string theory parameters 7 ~ (|5.20p ) and c is the loop 

production parameter. For ordinary four-dimensional strings c is extracted from simulations 
and is of order unity. For higher dimensional strings c is suppressed by a factor P given by, 

P~(-) (5.55) 

where 5 is the effective string thickness and d is the number of spatial dimensions. In our model 
5 ~ ry~^ ~ and d = Q giving, 

(5.50, 

and hence. 
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In higher dimensional theories the scahng parameter 7 (which determines the correlation length) 
is also expected to be suppressed by a factor of P during the scaling regime, but small scale 
structure on the strings is likely to lead to weaker P-dependence, P^s = f{P) [HIE]. However 
as the correlation length has already been determined directly in terms of the model parameters 
it is likely that this effect has already been accounted for. We also note that it appears consistent 
in the warped geometry to account for the warping of both space and time by introducing the 
transformations ti — )• ootj, r{ti) — )■ aQr{ti) into the expression ()5.53p above, yielding an additional 
factor of a^^. However it is our belief that this effect is already accounted for via the derivation 
of Ur in terms of the warped throat model parameters. In fact such an explanation gives a nice 
interpretation to the formula (j5.57p which is then seen as the product of three terms: a "warping 
term" which accounts for the back-reaction on the large dimensions a^"^, a term accounting for 
extra-dimensional effects {^/c^/R)^ and a standard Lorentz factor ~ g. 



The formation rate of black holes at t'{ti) is then equivalent to (minus) the rate of necklace 
formation at ti, 

dnpBH{t'{ti)) _ dn{r{ti),ti) ^ 3ur ,^ 
dt'{ti) ~ dti ^ ahf ' 

and finally we may use the the fact that 

dnpBH{t'{ti)) dnpBH{t'{ti)) dt'iU) 



dMpBH{t'{ti)) dt'{ti) dMpBH{t'{ti)) 

3l^r f dMT{ti)\'^ 



a2t. 



dti 



(5.59) 



to find the contribution to the PBH mass-spectrum from collapsing necklaces by substituting 
for t'{ti) and red-shifting to the current epoch. 



Of more interest for us is the calculation of the total contribution of the spectrum to the 
fraction of the critical density of the universe in the current epoch, ^IppHito)- The standard 
formula for ^pBH{to) from collapsing cosmic string loops is [20] 

npBH{to) = —^r dt^^^^^M{t,to) (5.60) 

Pcrit[to) Jmax{tc,U) "'^ 

where t^, is the formation time of a black hole with mass M^, « 4.4 x lO^'^/i"*^''^ gm (= lO^'' in 
Planck units), whose lifetime would be the present age of the universJ^. M(t, to) is the current 
mass of a black hole that formed at a time t, and tc is the time at which loops first begin to 
form. 



Again assuming that most of the loop production occurs in the radiation dominated era, the 
rate of black hole formation is then given by 

dnpBH ^ 3;^,./ a{tf , . 

dt ~ ah^a{toY 



^"^ Note that this value was calculated in a four dimensional FRW model using the standard value of G = 
6.67 X 10-"Afm2%-^ But in an a higher dimensional model, gravity is expected to become much stronger on 
very small scales resulting in a significantly higher rate of Hawking evaporation for the smallest PBH's. However 
we will neglect such small corrections. 
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where / is the fraction of loops which collapse on the first oscillation, which is expected to be 
small. It is typical to neglect the effect of Hawking radiation by making the approximation 
M(t, to) ~ M{t)^ which in the standard case {M[t) ~ oGt) has been shown to make a difference 
of less than 6 per cent to the final value of 0.pBH{to) LJ [55] . 

Adjusting the standard calculations to account for the effect of necklace formation, we expect 
there to be contributions to ^pBHito) from two qualitatively different sources: Firstly we expect 
to find a spike in the formation of PBH's in the very early universe due to the gravitational 
collapse of necklace loops which shrink to their minimum radius 5i ~ within a time t' ~ 
^1 + YQTx ) hough on cosmological time scales we may simply assume t' ~ tj. The total 

mass of all the beads contained in these loops is large, and is therefore the dominant contribution 
to the initial mass of the black hole. This corresponds to the first term in ()5.62p below. 

The second contribution to ^pBHito) comes from loops which collapse before shrinking to 
their minimum size, by adopting a sufficiently compact and spherically symmetric configuration 
on their first oscillation [TB]. In principle this process is continuous through the history of the 
universe, although at late times we may neglect the contribution of the beads to the masses of 
black holes formed in this way - leaving just the first term in the second integral. At early times 
the beads must be included and therefore both terms become important. 



The expression for the (approximate) contribution of PBH's formed from collapsing loops to 
the current mass-density of the universe is therefore 

npBHito) - f* du^^^^Mriu) X (Pi ' (Pi ' 

Pcrit{to) Jt+ \ ahf ahf 

where tf- are the times between which the necklace mass exceeds and tr marks the end of 
the radiation dominated era. In practice however the factor / ~ 10"^*' multiplying the integral 
on the second line indicates that (without fine tuning of the parameters) by far the largest 
contribution will come from the first integral between t~ and t^. In other words we expect the 
necklace-specific channel to dominate the production of PBH's and therefore choose to neglect 
the latter two terms. 



(5.62) 



The two factors in large brackets outside the first integral account for the red-shifting from 
the end of PBH production (with Mpbh > M*) from necklace collapse to the present day. For 
our purposes it is convenient to use the parameterisation 

/9erit(to) ~ Si^o'M^, ^ '^H^\ (5.63) 

Note that for PBH's to form via the necklace-specific process outlined above and to survive to 
the present day, thus contributing to the current mass density of the universe, we require 

Mriti) > (5.64) 

Though in our model the difference may be substantially greater and is therefore something which should be 
checked for completeness. 
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for at least some range of ti within tp^^ ^ ti < tp^^j. There are then three possible scenarios 
which we can consider. 

• < If this condition is satisfied for all ti in which black holes are produced, we can 
integrate the first term in (j5.62p over the range tpp^j < U < t~^pfj- 

• < Af* < With M^: in this range, the expression above will automatically be 
satisfied for times between < ti < t^, which may be obtained by solving (|5.64p . 

• > An PBH's formed by this process will evaporate long before the present 
epoch. 

In reality the first of these scenarios will not occur if 7]s ~ ^/a' ^ (= Mpi) unless R is hierarchi- 
cally larger than the fundamental string scale. This is theoretically compatible with the bound 
()5.12p for small enough values of the warp factor, though we need not assume that the bound 
is close to saturation. Working in Planck units and using the values for Q (j5.44p and rjs, it is 
possible to show that M^, > rj~^ /2Q so we need only 'tune' the values of our model parameters so 
that Mxiti) > for at least some range of ti. It turns out that the most important parameters 
for ensuring this condition is met are the warp factor oq and the world-sheet flux momentum 11. 
There is in fact a region (if somewhat restrictive) of the whole parameter space of ao,-R, 11 for 
which this condition is met which we shall discuss later. But it seems that one is either required 
to have small values of and/or large values of 11. That large 11 should help ensure PBH 
production is understandable since the pre-factor ^1 + simply multiplies the expression 

for MT{ti), whereas increasing world-sheet flux does not affect the value of M*. However the 
relation between Mxiti) and oq is more complicated as this factor appears in a complex way 
inside the Elliptic function. 

For the sake of completeness we calculate tf{M^) explicitly by setting Mr(tj) = M*. We 
then use the early time expansion of M^iti), given previously, to estimate t~ - and the late time 
approximation to Mxiti) to obtain t+. The result is 



4riaoa2ri^l + 5^ 



Although one could also obtain t+ by equating MTapprox{ti) with M^,, the resulting expression 
for t^ is rather unwieldy and does not significantly differ in its value from the approximate form 
given here. The evaluation of 0,pBH{to) (|5.62p between the limits t~-tf (given by ()5.65p above) 
is then obtained by integrating between t~ and tE using the early time approximation (|5.46p 
and between tE and tf using the numerical fit. The full expression is then well approximated 
by the following: 

Qito) ~ 
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-(2^^) ' - - -L(Log(256) - 3)(f2 - 1)1 + ^2^(1 + Log(16))(F2 - l)i 



1 



where 2-^1 is the usual Hypergeometric function and we have defined F = F{aQ, R, M^:,a, ATi, 11) 
to be; 

^ _ WalM^R^ 



As noted above, the fulfillment of the condition MT{ti) > M* requires small values of the warp 
factors and/or large values of the fiux parameter 11 - though the exact relationship between 
Mr(fj) and ao is complex. We must also ensure, for consistency, that t~{M^) < tf{M^) (where 
these values are given by the expression ()5.65p ). Exactly how small oq is required to be in order 
to fulfill both of these requirements depends on how large we take the 11 to be. We see from 
the explicit expression for 11 given below that, in principle, < IT < 00 for flux XFq^ in the 
allowable range < X'^F^^ < ag(a^r2 + R'^s''^). 

n = 2TiA^Fo^, 

al^r^ + a^'R^s'^-a^'X^Fl 

As causality requires that F^^ < a^A-2(r2 + a^^R'^s''^) we may assume, 

Fl = f34X-\r^ + a^^Rls'^) (5.68) 
where < /3 < 1. Equation (|5.67p may then be re-written in terms of (3 so that, 

from which it is clear that /3 = corresponds to 11 = and that 11 — t- 00 as /3 — )• 1. However 
the bound for PBH formation is not the only condition we must consider as the predictions of 
our model must also be consistent with observational bounds. Current observational constraints 
on the energy density of PBH's come from the EGRET experiment, which measures the extra- 
galactic gamma ray fiux at lOOMeV [27]. By calculating the expected contribution to this fiux 
from black holes expiring at the present epoch \38\ [39] (see also jlQ] for bounds derived for the 
standard PBH spectrum using the latest data) they were able to show that the current density 
of PBH's formed via the collapse of cosmic strings is bounded by 

i^PBHito) < 10"^ (5.70) 

Note however that this bound is based on the prediction that the PBH mass spectrum follows 
the profile predicted by the standard Hawking collapse process, such that MpBH{t) ~ t and 
the number density per mass interval is given by ()5.5ip . Technically one should recalculate this 
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bound using the spectrum predicted by the necklace-specific coHapse channel in order to place 
bounds on the model parameters from experimental data. We hope to supply a re-calculated 
bound on ^pBHito) ^ later date, but content ourselves for the time being with i^pBH{to) < 
10~^ as a "ball park" figure with which to proceed. 



With this in mind we see that not all the region of parameter space allowing PBH production 
is compatible with observation. For example, 'typical' values of oq ~ 10~^^ and IT ~ 10^^ are 
sufficient to ensure that Mr(tj) > M^, for at least some ti but the resulting value of the ^IpBHito) 
integral (|5.66p is huge {^IpBHito) ~ 10'^ >> 1) if R is comparable to the string scale. However 
due to the large iZ-dependence in the denominator caused by the extra-dimensional contribution 
to i^r we (happily) see that typical values of iZ ~ 10^ are sufficient to bring ^pBH{to) within 
the observable bound ^pBH{to) < 10~^. Thus the dimensions of the compact space need not 
be hierarchically larger than the string scale in order for the predictions of our model to be 
consistent with observational constraints. 



To systematically explore the values of aQ,R and flux 11 which are consistent with all the 
constraints above, we can express the value of the peak of the necklace mass function, Mxapprox 
given earlier in (|5.40p . in terms of M^,; 



n2 .,2^/1 



where k > 1 is a parameter that for given values of (oq, a, R, Yl, Ti, a') express how far the peak 
value of Mt is above the mass scale M^,. Using this definition we can expresses the quantity 



1 + in terms of (ao, a, R, k, Ti, a') and the above terms for tf simplify considerably 

0.171a'"^/^ , 2.297riKa'^/^ 



2rr ' * 2 



(5.72) 



It is interesting to note that both t~ and depend on the warp factor as Oq^. It may be 
thought that for k = 1, the times t~ and should approach one another since the definition 
of t^: requires values of ti for which Mt = . However recall we have assumed the early time 
approximation for Mt in determining t~ , and not the 'late-time' approximation MTapprox- This 
assumption then requires that we take k > 1. We could of course consider situations where k 
is closer to unity, but then one has to use MTapprox in the determination of both t~ and tf - 
which is technically complicated. 

Now evaluating (|5.66p using the above, we obtain an expression that depends on the param- 
eters {ao,a, R,Il,Ti,a' ,tQ, pcrit)- In Figure 8 we present a contour plot showing the values of 
R,K that are consistent with the condition ^psHito) < 10~^. In this plot we have taken, as 
a typical value, a = 0.8 and have set Ti = a'~^ = 1 and input the standard values of pcrit, to 
and M*. The resulting bounds on {R,k} are remarkably insensitive to the actual value of ao 
because the latter only appears in i^ppHito) via logarithmic factors. 

For given (allowed) values of {R, k} we can then deduce the value for the flux parameter 
n via (j5.7ip . which is of course sensitive to the value of oq. To visualize the allowed values of 
the flux parameter H corresponding to those of aQ,R,K - we use a 3D-contour plot shown in 



34 




35 



Figure 9. We illustrate six typical surfaces in the {Log[ao], R, k} plane corresponding to values 
of n = 10^^, 10^*^, 10^, 10^, 10^, 10^ (from front to back). These surfaces also take into account 
the observational allowed values of {R, k} illustrated in Figure 8. 

mo 
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Figure 9: 3D Contour plots showing values of the flux parameter 11 in the (Log[ao], R, k) 
plane. Six surfaces are shown corresponding (from front to rear) to the values of 11 = 
10^^, 10^'', 10^, 10^, 10^, 10^. These plots take into account the observational bounds on R,k 
shown in Figure 8. The remaining choice of parameters are as in Figure 8. 



5.3 Quasi-stable necklaces 

Let us now drop the assumption that the loops retain their necklace structure indefinitely. As 
outlined earlier, it is reasonable to expect that wrappings in the flat (/)-direction may unwind 
over time leading to a "flattening" of the ^-direction. This in turn flattens the ijj direction and, 
so on until the whole necklace structure unravels. 

Assuming that the motion of the string in the i;^)-direction remains random after the formation 
of a loop, and that the motion is also random in any newly flattened direction, we should expect 
the characteristic lifetime of the beads in any necklace to be comparable to the formation time, 
ti. This is because the average warped time taken for a single i;^- winding to contract to a point 
is ao^i ~ -? r-F - where this is the time associated with a single winding!^. The warped time 

^^Here we have used the fact that the radius of the winding must contract by a distance ~ R. Assuming 
a step length of ei ~ a^^\/a' for the change in the radial coordinate this requires a total displacement of 
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taken for n{ti) such windings to contract is therefore 



aotniu) ~ ^1 ^ (5.73) 



As we are only considering that part of the string which forms the extra-dimensional windings 
we may compare the expression above with (j5.3p in which ui and a have been set to unity. We 
then see that, 

tniu) ~ ti (5.74) 

This implies that the number of necklaces originally formed at ti which survive for a time 
At >> ti after formation will be negligible, with most having become standard string loops. 
However it is difficult to estimate the fraction of necklaces surviving for an arbitrary length 
of time At using such generic arguments, though this is exactly what we must calculate to 
determine the true contribution of necklace collapse to the PBH mass spectrum. In particular 
we must calculate the fraction of loops which retain their necklace structure for an interval 
At ~ pgy^ ti. To do this we consider the probability distribution which describes fluctuations 
of the radial coordinate. For a random walk this is simply a Gaussian distribution with mean 
H = and variance, 

a^{At) ~ flo ^^ef ~ a^'^Va'At (5.75) 

where St ~ a^^^/c^ is the time interval between steps, e; ~ Qq^Vc^ is the step length and At is 
the total (un-warped) time elapsed. The total probability density function ^{t, a{At)) is given 
below. 



"° xexpf--^). (5.76) 



x/2^(\/^At)2 V 2Va'At 

Since a loop forming at time tj has n(tj) windings to lose, the radial coordinate must "travel" a 
distance approximately equal to n{ti)R in order for the loop to lose its necklace structure. The 
fraction of necklaces /(At) which survive for an interval At after ti is then given by the integral 

of the above expression between t ~ ±aQ^n{ti)R ~ iag ^ [\/a'ti) ^ ~ ibcr(tj) which can be well 
approximated using the error function 




Thus when At ~ tj the fraction of loops which have retained their extra-dimensional windings 
is approximately 

fiU) ~ Erf (^-^) ~ 0.68 (5.78) 



.j^ steps. This in turn requires on average ~ /o-o ^ct' random steps which takes (a warped) time 



Clot ~ 3^ , X a,) = 
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which matches well with our result ()5.74p [fj. The fraction of necklaces which survive until they 
reach the minimum radius is, 

/ (At ~ -^t,) ^Erf{^[ — — ) ) (5.79) 




TGTi 

Strictly speaking therefore the integral (j5.62p should be multiplied by the additional numeric 
factor above to account for the the loss of necklaces which un-ravel before reaching the point at 
which they can undergo gravitational collapse. 

We also note that / (^At ~ j^*') ^ 10"^° when, 

a > 1.591 X W^^rgTi (5.80) 

and thus the necklace-specific channel becomes comparable (or sub-dominant) to the standard 
Hawking process. Practically however the values of a, T, Q and Ti here are such that this is 
unlikely to occur without significant fine tuning. 



5.4 Comparison with the standard string-monopole network 



The dynamics of string-monopole network evolution depends crucially upon the ratio of the 
monopole energy density to the string tension |17l I15j . This dimensionless parameter, denoted 
TZ in the literature is given by, 

7^=— (5.81) 
fid 

where m is the the monopole mass, /x is the tension and d is the distance between monopoles. 
For IZ « 1 the network behaves like an ordinary string network - reaching a scaling solution 
at late times, but for TZ» 1 the mass of the beads dominates the dynamics of the evolution. 
In our case these energy scales are time dependent and the corresponding parameter is given by 

= ^4r\ (5-82) 
^ ^ Tiditi) ^ ^ 

which scales as 7^ ~ for ti » and 7^ ~ t • ^ for ti « Thus we see that as 

tj — 7- oo we find 7^ — )• and vice-versa. Our results are therefore consistent with the standard 
analysis in that the necklaces behave like an ordinary string network for 7^ << 1 but like a 
string-monopole network when 7^ >> 1 since the mass of the beads becomes significant. 

However the dynamics of TZ differ profoundly in our model. This indicates that the effect of 
beads formed from extra-dimensional windings is different to the effect created by monopoles 
with regard to network evolution. To understand this in more detail let us consider the latter. 
The standard equation for the evolution of TZ is [41j . 

^ = -Kst-^ + Kgt-^ (5.83) 



^°Note that f{At) < 0.7 for all At > U. In fact Erf (1/2) ~ 1/2, hence the majority of loops (/(Af) > i) will 
have lost all their windings by At ~ 2ti. This helps to quantify our earlier result (|5.74|l more precisely. 
^^For example monopoles formed after a separate phase transition at some temperature Tm > Ts- 
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where the first term on the right hand side describes the stretching of the string due to the cosmic 
expansion and the the second term describes the contraction due to emission of gravitational 
radiation. It is a reasonable to assume that 



K = Kg - Kg > (5.84) 
which allows us to solve (j5.83p . up to some constant of integration, 

7^ ~ t". (5.85) 
Using (j5.8ip the evolution of the inter-monopole distance then scales like 

d~7^-l~^-'^. (5.86) 

This was the crucial assumption Matsuda used to identify d{t) with the step length of the random 
walk x(t) |41j of the monopole along the string. From this point, the estimated initial number 
of beads per loop is then given by 

^-^^^^n-i (5-87) 



d{ts) X 

where d{ts) ~ (tjvf^s)^ is the initial bead spacing in the network at the time of string formation 
tg (where here tg ~ Vc/) and k = corresponds to the (natural) k = 1 solution of ()5.83p . 

Assuming that ~ ^/n beads would survive without annihilation, we find that nicoii oc ^/n 
and hence, 

rucoii ~ tf^ (5.88) 



which favours black hole production at late times vs dark matter production for small ti. This 
is the standard result assuming that necklaces which reach their minimum radius but which 
have insufficient mass to undergo collapse to form PBH's can only interact with other matter 
gravitationally, leading to the production of Planck scale dark matter relics. In Matsuda's 
original scenario however, based on the assumptions above, large numbers of Planck scale dark 
matter 'particles' are formed at early times with smaller numbers of increasingly massive black 
holes formed at later times. In our model this process is essentially reversed with a window 
of PBH production in the early universe (though potentially this window may be quite large 
for small oq so the last epoch of PBH formation is not necessarily 'early', even on cosmological 
time scales) when Mxiti) > with dark matter forming well into the scaling regime when 



Physically the above equations describe the shrinking of the string sections connecting neigh- 
bouring monopoles. Effectively the contraction of the string is able to pull these monopoles 
through the horizon at an ever increasing rate El. Hence we find that the early and late time 
limits of TZ as defined above, are the opposite of those obtained from ()5.82p . namely that 7?. — )• oo 
as time increases, whilst tending to zero as ti — )• 0. 

^^Alternatively this can be viewed as the expanding horizon uncovering monopoles, separated by increasingly 
short distances, the distances having been shortened by the contraction of the string. 
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Similarly the natural solution, k = Kg — Kg ^ 1 obtained from an order of magnitude 
estimate, gives d ^ in the standard case as opposed to d ~ t in our case. In yet another 
respect therefore we have obtained the opposite results to the standard analysis. 

However it is not immediately obvious that what applies to monopoles formed in a sepa- 
rate phase transition also applies to beads formed by extra-dimensional windings. Although 
a contracting string may pull ordinary monopoles through the horizon at an ever increasing 
rate, a winding may not cross the horizon ready-made and concentrated at a point from a 
four-dimensional perspective. This is because insufficient time will have elapsed to establish 
correlations in the compact dimensions. Put another way, if the horizon advances by a distance 
aoc5t, the end point of the string formerly localised at the horizon, can move by, at most, a 
distance c6t from its original position in the compact space. Causality therefore places a limit 
on the rate at which new windings can enter the horizon, whereas no such limit exist for ordinary 
monopoles El. 

In fact the considerations above help to explain why the number of windings n(tj) is pro- 
portional to the ratio ao/R not simply to in the warped throat model. Although larger 
R results in a slower rate of winding formation (as we would expect) this effect is, at least 
potentially, dwarfed by the effect of the warp factor (c.f. ()5.17p ) which is related to R via 
the deformation parameter, ~ e^R~^ This is because ao limits the increase of the horizon 
distance in the infinite directions (but not in the compact space) via 

dH{5t) = aoc5t (5.89) 

As the resultant velocity of the end point of the string is c = 1 smaller ao results in a greater 
velocity in the compact spaces and hence a larger limiting value for n{ti). 

In general however, windings may enter the horizon more rapidly the smaller their effective 
radii. In fact this at least partially accounts for the falling bead mass at late times; since 
tj — 7- oo we expect that tj/ — t- for a random walk in the extra dimensions. This could be 
achieved through either; 

• Falling number density (per unit distance along the string) at late times - n ~ 

• Windings wrapping ever smaller effective radii in the 



The fact that d ~ tj at late times highlights the first effect and M;, ~ t^ indicates the 
second. In practice however we would expect both effects to play some role. 

The question then arises, does the shrinking bead mass have a counterpart in field theory? 
We may also consider whether the arguments for the limitation of windings entering the horizon 
hold true for field-theoretic strings in the dual non-Abelian Higgs models. The present authors 
hope to address this question in a future publication. 

^^An alternative way to view this result is to consider windings as correlations in the compact space. All 
windings must therefore form within the horizon to preserve the causal structure. 
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Finally we note that the predictions of our model (specifically that PBH production is 
favoured at early times with dark matter production at later epochs, in contrast to previous 
predictions) are highly dependent on the assumption of a random walk regime which produces 
the extra-dimensional windings. It is not immediately clear how this behaviour might be modi- 
fied by the introduction of a velocity correlations regime, though it is possible that the results 
may be more in line with those of previous studies. This is a question the authors hope to 
pursue at a later date. It is also highly dependent upon the assumption that only beads formed 
from net windings, not from bead-anti-bead (6 — b) pair formation, contribute to the mass of 
PBH's/DM relics. As we have seen, assuming that ~ ^/n 6/6's survive until the string reaches 
it's minimum radius (where n here must be taken to include all beads/anti-beads created from 
random steps along the string, not simply those formed by net displacement in the compact 
space) leads to ()5.88|) and to qualitative behaviour which is the very opposite to that obtained 
in our model. However as stated previously in the footnote on page 15, we believe that such an 
argument is valid only for a static string and that in a contracting loop all 5's will eventually 
collide with their b counterparts before the minimum radius is reached. In particular, in our 
model we may expect the initial b — b spacing to be approximately equal to the step length 
ei ~ ttQ^ay/c^ which is most certainly greater that the minimum radius 6i ~ r]~^ ~ y/o/ for all 
reasonable values of a and ao- 



6 Discussion 



In this note we have investigated a simple cosmic string model using several key ideas from 
string theory. Our starting point was the assumption that there is some initial stage of inflation 
that leads to the creation of trapped strings (along the lines of [7]) that have non-zero windings 
in the internal space. We argued that the string feels a lifting potential due to the non-trivial 
geometry of the background, forming a 'kinky cycloop'. By estimating the mass of such a string 
configuration we were able to show that formation at late times reduces to a scaling solution, 
which one may hope to identify with a relic dark matter phase. We also note in passing that, in 
principle, we require dark matter production from necklaces between f^^jj and to not to over- 
close the universe. However, as the present day dark matter density bound is approximately 
^Dhlito) < 0.3 and due to the t~^ suppression of the loop production function, we find that in 
all but pathological cases the PBH density bound Op_B//(to) < 10^^ is by far the most stringent 
constraint. 

In some sense this has been a bottom up approach to the problem. The warp factor is simply 
a parameter of the metric, depending only on the extra-dimensional flux and the deformation 
parameter of the (non-compact) supergravity background. A fully realised UV approach will 
precisely fix the warping as a function of the closed string moduli VEV's and the flux, therefore 
leading to a more tightly constrained theory. However we would hope that the general features 
of our approach will remain intact. In particular the formation window for PBH at early times 
contrasts strongly with the field theory/string model proposed by Matsuda. Although the 
authors initially feared that result may be in some sense disappointing, due to the likelihood of 
PBH's formed at early times having evaporated by the present epoch, we have shown that it is 
possible within the model here constructed to satisfy the condition MT{ti) > and current 
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observational constraints. It is worth noting too at this point that, for very smah oq, the window 
of black hole formation may be sufficiently large for the assumption of a constant M* to become 
inaccurate. A more detailed analysis would therefore take into account the rate of black hole 
evaporation in the extra-dimensional background introducing a time-dependent lower bound 
on the PBH mass M*(tj). An improved analysis would also consider how the abundance of 
PBH's formed through necklace collapse in the very early universe affects BBN and subsequent 
structure formation. In fact, there are many theoretical and observational constraints which 
any such model must satisfy concurrently, thereby potentially allowing the Klebanov-Strassler 
model parameters to be bounded with even greater precision. 

It is clear from our model that in order to generate any PBH's of mass M* restricts us to 
a small region of the {ao,R,Il) parameter space as is shown by in Figures 8 and 9. This is not 
surprising if we recall that M^, ~ 10'^'' in Planck units and so there is a hierarchy issue here. On 
the other hand since the warp factor oq depends exponentially on background fluxes, obtaining 
small enough values to generate a hierarchically large Mt is perhaps not so unnatural. 

In order to satisfy the PBH observational bounds one requires larger background flux M if 
there is larger world- volume flux H. This is inherently obvious, since the additional world- volume 
flux can be treated as an effective mass correction to the necklace. This additional mass will 
back-react on the solution and the probe limit will be rendered invalid, unless the background 
flux is also increased to compensate. What is interesting is that, for fixed warping, there is a 
large parameter space where bounds are satisfied. Another interesting feature is that, as is clear 
from Figure 8, the allowed values of R for regions where Mj- >> (corresponding to large k) 
are still in the 'mild' hierarchy region R ~ 10^/^. Indeed the only way to have a value of i? ~ 
and consistent with observation requires a fine tuning where the peak of Mt is taken very close 
to but just above M^,. This then reduces the time interval of mass M^, PBH formation given by 
(t^ ~ tt) which can be made small enough so that we can reduce the value of R and still obtain 

npBHito) < w-^- 

Our construction is very model dependent in some sense, since it is only valid for a special 
class of supergravity backgrounds of the type IIB string. Another interesting class of cosmic 
string models arises from considering Heterotic M-theory on the orbifold CY3 x /Z2- The 
strings in this case can arise from wrapping membranes (or five-branes) over various cycles within 
the internal space |42tl43|. Topologically stable strings are only possible for five-branes wrapping 
a complete four-cycle within the CY space [S]. However, other wrappings are potentially 
possible if there is an uplifting potential and a similar analysis can be performed. We leave this 
for a future publication. 

One other issue that remains is the validity of substituting the expression r(tj) = ati for the 
variable r in the bead mass terms of the EllipticE expansion and r{t,ti) = ati — TTiQ{t — ti) 
into the leading order term representing the mass of the four-dimensional sections of the string. 
As stated previously, this is entirely consistent with Matsuda's original assumption that the 
bead mass remains constant after loop production due to the "trapping" effects of the lifting 
potential, whereas the sections of string not contained in the extra-dimensional windings shrink 
through the emission of gravitational radiation. However, the re-paramaterisation invariance 
of the Nambu-Goto action implies that it may not be possible to talk of one 'section' of the 
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string expanding/contracting while other sections remain fixed. By definition strings have no 
sub-structure which would seem to imply that they must expand or contract along their entire 
length. In our model this would mean that the contraction of the string due to gravitational 
radiation in the warped Minkowski directions in fact produces a contraction of the string along 
the internal direction as well. If true this implies that the expression r{t,ti) should be substituted 
in both the bead-mass and four-dimensional parts of the string mass integral, thus creating an 
explicitly time dependent bead-mass Mb{t,ti) for t > ti. Potentially this may alter the present 
results significantly, as we would expect the number of necklaces which retain sufficient mass 
in their extra-dimensional windings to produce PBH's when reaching their minimum radius 
to diminish sharply. This is therefore likely to lead to a much smaller amount of black hole 
production and greatly increased dark matter formation. It is hoped that the question of the 
explicit time dependence of the bead mass may be resolved conclusively by investigating field 
theory duals of the string necklaces discussed here. Again, the present authors hope to be able 
to offer some contribution in this direction in the near future. 

A further aside concerns the question of binding energy between b — b oi b — b pairs which 
coalesce. In field theoretic models the binding energy between true monopole or anti-monopole 
pairs may be calculated. We may suspect that, if a similar binding energy exists between beads 
formed from extra-dimensional windings, it's presence may alter the string dynamics so that 
the evolution of a string configuration with n 'one step' beads into a single bead formed from 
n steps will be energetically favoured. Once formed it is also possible that such a configuration 
may un-wind more rapidly than it's corresponding n-head counterpart, leading to a necessary 
modification of the arguments in Section 5.3. As a first approximation we have assumed that 
there is no binding energy between winding states. In the absence of a full BSFT picture this 
seems intuitively likely, given that the string must be considered as a collective phenomenon, 
which suggests that the unwinding of n single steps will be just as hard as un-winding a single 
n-step wrapping. However it is not clear how one can really understand this problem without a 
more detailed analysis of the boundary theory. 

Finally it is interesting to speculate further on the recent result implied by the PAMELA 
experiment |451 H6] which suggests that there is an over-abundance of high energy positrons in 
the cosmic background. Rather than interpreting this as a signal for dark matter, it has also 
been suggested that this is in agreement with predictions from cosmic-strings |47] . It would 
certainly be interesting to extend our analysis along those lines. 
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